Astrometric and Timing Effects of Gravitational Waves from Localized Sources 



Sergei M. Kopeikin* and Gerhard Schafer 
Theoretisch-Physikalisches Institut der Friedrich-Schiller-Umversitat Jena, Max-Wien-Platz 1, 077^3 Jena, Germany 

Carl R. Gwinn 

Department of Physics, Broida Hall, University of California at Santa Barbara, Santa Barbara, CA 93106 

T. Marshall Eubanks 
US Naval Observatory, 3450 Massachusetts Avenue, Washington, DC 20392 



Contents 



0^ 



> 
m 
o 
o 



00 = 
On . 

I 



X 



Introduction 



II 



A Historical Rcmarkf 



B Observational Capabilities 



Equations of propagation of electromagnetic waves 



III Metric Tensor and Coordinate Systems 



IV 


Method of Integration of the Equations of Motion 








A 


Useful Relationships 








B 


Calculation of Integrals from the Static Part of the Gravitational Field] . 






C 


Calculation of Integrals from Time Dependent Part of Gravitational Field 






1 First Method of Integration 




2 Second Method of Integration! 



3 

3 

4 

4 
5 

8 

8 



V Perturbations of Photon's Trajectory 



15 



VI Basic Observable Relativistic Effectg| 17 

[A Time Delay]. . . ' 18 

B Deflection of Light] 19 



VII Discussion 




A Case 1. Small Impact Parameter frn < 0) 




1 Asymptotic expansions of independent variables 


2 Asymptotic expressions for time delay and the angle of light deflection 




B Case 2. Small Impact Parameter (rn > 0) 




1 Asymptotic expansions of independent variables 


2 Asymptotic expressions for time delay and the angle of light deflection 




C Case 3. Large Impact Parameter 




1 Asymptotic expansions of independent variables 


2 Asymptotic expressions for time delay and the angle of light deflection 





22 



VIII Conclusions 



APPI NDIXES 



A Comparison to the paper by Damour &; Esposito-Farese 



B Harmonic and ADM gauge conditions in the first post-Minkowskian approximation 



29 
30 
30 



33 
33 
33 
35 



*On leave from: ASC FIAN, Leninskii Prospect, 53, Moscow, 117924, Russia 



1 



The extremely high precision of current radio interferometric observations demands a better 
theoretical treatment of secondary effects in the propagation of electromagnetic signals in variable 
gravitational fields. Such fields include those of oscillating and preccssing stars, stationary or coa- 
lescing binary systems, and colliding galaxies. Especially important is the problem of propagation of 
light rays in the field of gravitational waves emitted by a localized source of gravitational radiation. 
A consistent approach for a complete and exhaustive solution of this problem is developed in the 
present paper in the first post-Minkowskian and quadrupole approximation of General Relativity. 
This approximation is linear with respect to the universal gravitational constant G and accounts 
for the static monopole, spin, and time-dependent quadrupole moments of an isolated system. We 
demonstrate for the first time that the equations of light propagation in the retarded gravitational 
field of an arbitrary localized source emitting quadrupolar gravitational waves can be integrated 
exactly in closed form. The infiuence of the gravitational field under consideration on the light 
propagation is examined not only in the wave zone but also in cases when light passes through 
the intermediate and near zones of the source. We reproduce the known results of integration of 
equations of light rays, both in a stationary gravitational field and in the field of plane gravitational 
waves, establishing the relationship between our new formalism and the simplified approaches of 
other authors. Explicit analytic expressions for light deflection and integrated time delay (Shapiro 
effect) are obtained accounting for all possible retardation effects and arbitrary relative locations of 
the source of gravitational waves, the source of light rays , and the observer. Coordinate dependent 
terms in the expressions for observable quantities are singled out and used for physically meaningful 
interpretation of observable quantities. It is shown that the ADM and harmonic gauge conditions 
can both be satisfied simultaneously outside the source of gravitational waves. Such ADM-harmonic 
coordinates are extensively used in the present paper. Their use drastically simplifies the integra- 
tion of light propagation equations and the equations for the motion of light source and observer 
in the gravitational field of the source of gravitational waves, leading to the unique interpretation 
of observable effects. The two limiting cases of small and large values of impact parameter d are 
elaborated in more detail. It is proved that leading order terms for the effect of light deflection in 
the case of small impact parameter depend neither on the radiative part (~ of the gravitational 
field nor on the intermediate (~ zone terms, confirming a previous result in the literature. The 

main effect rather comes from the near zone terms. This property of strong suppression of 

the influence of gravitational waves on the propagation of light rays makes much more difficult any 
direct detection of gravitational waves by VLBI or pulsar timing techniques, in contrast to previous 
claims by other authors. We also present a thorough-going analytical treatment of time delay and 
bending of light in the case of large impact parameter. This exploration essentially extends previous 
results regarding propagation of light rays in the field of a plane monochromatic gravitational wave. 
Explicit expressions for Shapiro effect and deflection angle are obtained in terms of the transverse- 
traceless (TT) part of the space-space components of the metric tensor. We also discuss the relevance 
of the developed formalism for interpretation of radio interferometric and timing observations, as 
well as for data processing algorithms for future gravitational wave detectors. 
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I. INTRODUCTION 



A. Historical Remarks 

Binary systems are well known sources of periodic gravitational waves. Indirect proof of the existence of gravitational 
waves emitted by binary pulsars was given by Taylor However, the direct observation of gravitational waves still 
remains the unsolved problem of experimental gravitational physics. The expected spectrum of gravitational waves 
extends from - 10'' Hz to lO-^^ Hz |], |]. Within that range, the spectrum of periodic waves from known binary 
systems extends from about 10~^ Hz, the frequency of gravitational radiation from a contact white-dwarf binary [Q, 
through the lO""* to 10~^ Hz range of radiation from main-sequence binaries [|||, to the 10~^ to 10~^ Hz frequencies 
emitted by binary supermassive black holes postulated to lie in galactic nuclei^] . The dimensionless strain of these 
waves at the Earth, h, may be as great as 10^^^ at the highest frequencies, and as great as 3 x 10"^^ at the lowest 
frequencies in this range. 

Sazhin ^ first suggested detection of gravitational waves from a binary system using timing observations of a pulsar, 
the line of sight to which passes near the binary. It was shown that the integrated time delay for propagation of the 
electromagnetic pulse near the binary is proportional to where d is the impact parameter of the unperturbed 
trajectory of the signal. More recently, Sazhin & Saphonova made estimates of the probability of observations of 
this effect, for pulsars in globular clusters, and showed that the probability can be high, reaching 0.97 for one cluster. 
We note however that mathematical technique worked out in these papers allows rigorous treatment only of effects 
of the near-zone, quasi-static quadrupolar part of the gravitational field and is not enough to make any conclusion 
about actual observability of gravitational waves emitted by a binary system. 

Wahlquist |^ made another approach to the detection of periodic gravitational waves, based on Doppler tracking 
of spacecraft traveling in deep space. His approach is restricted by the plane gravitational wave approximation 
developed earlier by Estabrook & Wahlquist jl^ . Tinto ( [0 , and references therein) made the most recent theoretical 
contribution in this area. The Doppler tracking technique has been used in space missions, by seeking the characteristic 
triple signature, the presence of which would reveal the influence of a gravitational wave crossing the line of sight 
from spacecraft to observer JT^ . 

Quite recently, Braginsky et al. [|l3| (see also [Q) have raised the question of using astromctry as a detector of 
stochastic gravitational waves. This idea has also been investigated by Kaiser & Jaffc and, in particular, by Pyne 
et al. and Gwinn et al. |l7| who showed that the overall effect is proportional to the strain of metric perturbation 
caused by the plane gravitational wave and set observational limits on the energy density of ultra long gravitational 
waves present in early universe. Montanari studied polarization perturbations of free electromagnetic radiation 
in the field of a plane gravitational wave and found that the effects are exceedingly small. 

Fakir ( |^ , and references therein) has suggested the possibility of using astrometry to detect periodic variations in 
apparent angular separations of appropriate light sources, caused by gravitational waves emitted from isolated sources 
of gravitational radiation. He was not able to develop a self-consistent approach to tackle the problem with necessary 
completeness and rigor. For this reason his estimate of the effect is erroneous. Another attempt to work out a more 
consistent approach to the calculation of the deflection angle in the field of arbitrary source of gravitational waves has 
been undertaken by Durrer [po[ . However, the calculations have been done only for the plane wave approximation and 
the result obtained was extrapolated for the case of the localized source of gravitational waves without justification. 
For this reason the deflection angle was overestimated. The same misinterpretation of the effect can be found in the 
paper by Labeyrie pH who studied a photometric modulation of background sources by gravitational waves emitted 
by fast binary stars. Because of this, the expected detection of the gravitational waves from the observations of the 
radio source GPS QSO 2022-1-171 suggested by Pogrebenko et al. was not based on firm theoretical ground. 

Damour & Esposito-Faresc [|2^ have studied the defiection of light and integrated time delay caused by the time- 
dependent gravitational field generated by a localized material source lying close to the line of sight. They explicitly 
took into account the full, retarded gravitational field in the near, intermediate, and wave zones. Contrary to the 
claims of Fakir ||l^ and Durrer |^ and in agreement with Sazhin's calculations, they found that the deflections 
due to both the wave-zone gravitational wave and the intermediate-zone retarded flelds vanish exactly. The leading 
total time-dependent deflection is given only by the quasi-static, near-zone quadrupolar piece of the gravitational 
field. 

In the present paper we work out an even more systematic approach to the problem. While Damour & Esposito- 
Farese considered both the light source and the observer to be located at infinity, and performed their calculations 
in terms of the spacetime Fourier transform, we do not need these assumptions. Our approach is much more general 
and applicable for any location of the source of light and observer in space with respect to the source of gravitational 
radiation. The integration technique which we use for finding the solution of equations of propagation of light rays 
was partially employed in [M and does not require any implementation of the spacetime Fourier transform. 
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Section 2 of the present paper discusses equations of propagation of electromagnetic waves in the geometric optics 
approximation. The metric tensor and coordinate systems involved in our calculations are described in section 3 
along with gauge conditions imposed on the metric tensor. The method of integration of the equations of motion with 
emphasis on specific details of calculations of particular integrals is given in section 4. Exact solution of the equations 
of light propagation and the form of relativistic perturbations of the light trajectory are obtained in section 5. Section 
6 is devoted to derivation of basic observable relativistic effects - the integrated time delay and the deflection angle. 
We find the more precise meaning of quite general formulae obtained in the previous section by discussing in section 7 
several limiting cases in the relative configuration of the source of light, the observer, and the source of gravitational 
waves. Section 8 contains concluding remarks. Appendix A compares results of our calculations with those by Damour 
& Esposito-Farese and proves their gauge invariance. Appendix B gives more details on the derivation of the 
ADM- harmonic coordinate system used in the present paper for interpretation of observed relativistic effects. 



B. Observational Capabilities 



Calculations of the effects of gravitational waves are of most interest if they indicate that those can be detected with 
present techniques, or foreseeable improvements. Astrometric precision and accuracy have evolved rapidly in the last 
decades, and can be expected to continue to improve. In principle, the accuracy attainable with a given instrument 
is approximately the angular resolution of the instrument, divided by the signal-to-noise ratio. In practice, the time 
span of the observations and the angular separation of the source from reference sources critically affect the attainable 
accuracy. 

Very-long baseline interferometry, or VLBI, attains the highest angular resolution available on an operational basis. 
It achieves angular resolution set by the diffraction limit, of « A/-B, where B is the separation of the interferometer 
elements (the baseline), and A is the observing wavelength. Practical baselines may be about as long as an Earth 
radius, B ~ 6400 km; a typical observing wavelength is A = 3 cm, yielding angular resolution of 1 milliarcsecond. 

Observations of a moderately strong (~ 1 Jy) extragalactic source, such as a quasar, can reach signal-to-noise ratio 
of several hundred in 5 or 10 minutes, offering potential angular accuracy of microarcseconds. In principle, a day of 
integration with the US Very Long Baseline Array (VLBA) could yield angular accuracy of about 0.1 microarcseconds. 

Observations using the largest radiotelescopes can increase the signal-to-noise ratio by a factor of ~ 10 x. In 
practice, a host of geodetic and propagation effects limit the reproducibility of VLBI astrometry. These factors must 
either be measured during the observations, or calculated from models. At present, atmospheric stability and changes 
in source structure limit reproducibility of measured angles between sources to about 1 milliarcsecond, over periods 
of months. Observations of pairs of radio sources, with separations of 0.5°, can yield angular accuracy of about 
50 microarcseconds, reproducible over periods of years, when effects of source structure are included 

Astrophysical H2O masers have extremely high flux densities, of up to 10^ Jy at A = 1.3 cm. In principle, a day of 
observation of masers with the VLBA could yield angular accuracy of a few picoarcseconds. Observations of masers 
have attained reproducibility of better than 10 microarcseconds over several months, between individual maser spots in 
a Galactic maser cluster, with separations of a few arcsec ||2^. Astrometric observations of extragalactic masers have 
attained accuracies of better than 1 microarcsecond, for maser spots separated by less than 1 arcsec 1 27 . Atmospheric 
variations probably dominate the error budget. 

Shorter wavelengths offer potentially higher diffraction-limited angular resolution, but practical obstacles are severe. 
Atmospheric effects present greater phase changes, on shorter time scales; and photon shot noise becomes a limiting 
factor for fainter sources and at shorter wavelengths. Optical interferometers in space will probably equal and exceed 
the accuracy of VLBI. For example, the Space Interferometry Mission (SIM), and the proposed European mission 
GAIA seek to attain angular accuracy of about 1 microarcsecond in several hours of integration ||2^ - Q . 

Astrometric observations to seek effects of gravitational waves could attain higher accuracy, at least on shorter 
timescales. The periods of the waves, and of the expected deflection, could be short enough to avoid some atmospheric 
and other propagation effects. For known binary systems, the wave period, and perhaps the phase, are known 
accurately, permitting search for deflections at this period. Such a "synchronous" search would eliminate many noise 
sources, allow detection of short-period motions with the sensitivity resulting from long integrations, and perhaps 
allow astrometric accuracy to approach the signal-to- noise ratio limit. 



II. EQUATIONS OF PROPAGATION OF ELECTROMAGNETIC WAVES 



We assume the approximation of geometric optics, as the wavelength of electromagnetic waves used for astrometric 
observations is usually much smaller than wavelength of gravitational waves emitted by isolated astronomical systems 
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like binary stars or supernova explosions This allows us to use the relativistic equation of geodesic motion of a 
massless particle (such as a photon) for description of the process of propagation of electromagnetic signal from the 
source of light to the observer at the Earth. We also assume that space-time is asymptotically flat. This assumption 
does not hold for cosmological distances. However, if we neglect all terms depending on the rate of cosmological 
expansion and make a rescaling of time and space coordinates with the cosmological scale factor a(t), our formalism 
will be still valid for application in cosmology. 

We denote spatial coordinates by = x = {x^,x'^,x^) and time coordinate — ct, where c is the speed of light 
and t is coordinate time. Let the motion of a photon be defined by fixing the mixed initial-boundary conditions 
introduced and extensively used by Brumberg 

/ N (ix(— oo) , 

X to = xo, ^ ' = k, 1 

at 

where = 1 and the spatial components of vectors are denoted by bold letters. These conditions define the 
coordinates xq of the photon at the moment of emission Iq and its velocity at the infinite past and the infinite 
distance from the origin of the spatial coordinates (that is, at past null infinity). In what follows we put c — I for 
convenience. 

Equation of propagation of photons in a weak gravitational field is given in the first post-Minkowskian approximation 
by the formula |3|,^: 

^'(t) = ^9oo,i - 9oi,t - ^goo,ti' - gik,ti^ - {goi,k - 9ok,i) i'' - (2) 

gmM^^x" - {^ik,] - ]^gkj,i^ x^x^ + Qff/ci,t - ffofej^ x''x^i\ 

where the goo, goi, dij are components of metric tensor, fully determined by the given distribution and motion of 
mass inside the source of gravitational field, dots over vectors denote the total derivative with respect to time, and 
commas indicate partial derivatives with respect to spatial or time coordinates; that is, for any function f^i = df /dx^, 
f^t = df /dt. Hereafter repeated latin indices mean summation from 1 to 3. The given equation is valid in arbitrary 
coordinates (gauges) and represents the ordinary second order differential equation for light propagation. 

The right-hand side of equation (^ includes terms which depend on the coordinate velocity of the photon, in 
the weak-field approximation approximately equal to the speed of light c. We restrict ourselves to finding a solution 
of equation (^ only in the first linear approximation with respect to the universal gravitational constant G. For this 
reason, when solving equation (^), only one iteration is enough and it is admissible to make the replacement = fc' 
in the right-hand side of the equation. The result of this approach is: 

x'it) = ^goo,i - 9oi,t - ^5oo,tfc* - gij,tk^ - {gozj - goj,i) - (3) 



ioo.j^'' (^gipj 2^pj;i^ + (^2^pj;t ~ gop.j^ k^k-'k^ 



This equation must be solved to obtain a perturbed trajectory of the photon propagating through the gravitational field 
of an isolated astronomical system emitting gravitational waves. To accomplish this task one needs a mathematical 
expression for the metric tensor. 



III. METRIC TENSOR AND COORDINATE SYSTEMS 



Let us chose the origin of the asymptotically fiat coordinate frame at the center of mass-energy of the isolated 
astronomical system and impose the de-Donder (harmonic) gauge conditions on components of the "canonical" metric 
tensor. We assume that gravitational field is weak and the metric of spacetime gap is written as a sum of the Minkowski 
metric rjajs = diag(— 1, 1, 1, 1) plus a small perturbation /iq,^: 



gal3 = r]af3 + hafl, (4) 

where the Greek indices run from to 3. The most general expression for the linearized metric tensor, generated 
by a system emitting gravitational waves, in terms of its symmetric and trace-free (STF) mass and spin multipole 
moments is given by Thorne Q (see also [^,|5|). It can be written as 
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ha0 = h"^^- + V^Wa + VaW^ , (5) 

where — The "canonical" form of the metric tensor perturbations in harmonic gauge reads as follows js^ 



'00 



(-1)' 



lean. 

''■Oi — 



2M 

— + 2E , 



(6) 
(7) 



■LCa7l. 



1=2 



(2M 



^ipqSpAi-l {t — r) 



1=2 



,qAi-i 



^(-ly [JA,(i-r) 



4E 



1=2 



n 



(8) 



^(^ + 1)! 



where TV* = a;*/'', and the round brackets around indices in equation (^) means symmetrization; that is, for any 
two indices T(^ij) = ^{Tij + Tji). In the pure harmonic gauge the functions , are solutions of the homogeneous 
d'Alembert's equation and are given by the expressions 



1=0 



(9) 



E 

/=o 



,iAi 



E 

1=1 



ZqAi^At- 



(10) 



where Wa, , '^Aj, 3^iA;_i, and Z^Ai^i are arbitrary functions of time. Their specific choice will be made later on 
in the discussion regarding the interpretation of observable effects. In equations (|^)-(10), we adopt the notation: 
Ai = aia2-.-ai is a polyindex, Ai is the total (Tolman or ADM) mass of the system, Tai and Sai are the STF mass 
and spin gravitational multipoles, and Wa, , Xai , 3^a, , 2a, are multipoles which reflect the freedom of coordinate 
transformations. These multipoles can be eliminated from the metric using the transformation 

a;'" = - , (11) 

relating an original harmonic coordinate system to another harmonic one x'", in which only the "canonical" part 
of the metric is present. 

However, we would like to emphasize that, in general, equation holds in an arbitrary gauge. Particular examples 
of functions w'^ and in harmonic gauge are given in equations (^-(^|). Other expressions for and in the 
ADM (Arnowitt-Deser-Misner) gauge are given in Appendix B wherein we also prove that it is possible to choose 
functions w*^ and in such a way that ADM and harmonic gauge conditions will be satisfied simultaneously. This 
means that the classes of harmonic and ADM coordinates overlap. The discussion of different gauges is helpful for giv- 
ing a unique interpretation of observable effects by properly fixing the coordinate degrees of freedom in corresponding 
quantities p8[ . 

The STF cartesian tensor has a special algebraic structure which eliminates all reducible parts of the tensor and 
leaves only the irreducible part having the highest rank |^,^. In other words, contraction over of any two indices of 
STF tensor gives identically zero. It is worth noting the absence of the dipole mass multipole 2i in equations (^)-(||) 
which is identically zero, due to the choice of the origin of coordinate system at the center of mass of the gravitating 
system. We also stress that the multipoles in the linearized metric (^)-(lO) depend on the "retarded time" t — r. At 
first sight this dependence seems to make subsequent calculations more difficult. However, just the opposite happens 
and the dependence of the multipoles on the retarded time makes the calculations simpler. 



6 



In what follows we consider the concrete case of a localized deflector emitting gravitational waves. In this section 
we restrict ourselves to considering the influence of gravitational field of the deflector on the propagation of electro- 
magnetic signals made by its total constant mass M, spin S, and time-dependent quadrupole moment Tij {t — r) only. 
This simplifies the expressions (^-(^ for the metric tensor, which are reduced to the expressions 



i^can. 

"-00 — 



2M 



Ipq(t-r) 



(12) 



i^can. 



2V, 



(13) 



T^can. f lean. 

'Hj — "ij "'00 



(14) 



where 



2 ■■ 

(t - r) 



(15) 



Herein terms depending on and Si are static and produce well-known effects in the propagation of light rays. 
Retarded terms that are inversely proportional to the distance r from the gravitating system describe the pure 
gravitational-wave part of the metric. 

Let us stress that in the harmonic coordinate system the gravitational-wave part of the metric tensor is present in 
all of its components and is expressed through the second time derivative of the quadrupole moment [ p3[ . If we choose 
the ADM gauge it is possible to eliminate the gravitational wave terms from the /ifjo"' and /i^""" components of the 
metric tensor and to bring all of them to h""^"^ | . Then /iqo™ ^''^^ ^of"^ depend only on the "instantaneous" time t 
and not on the retarded time t — r (see Appendix B). In combining the ADM gauge with the harmonic gauge an even 
simpler representation is possible where hoo and hoi do not depend on time at all. However, the transformation from 
the canonical form of metric ([T^)-([r^) to the ADM- harmonic form includes integration of the quadrupole moment 
with respect to time. Appendix B gives a more detailed study of this procedure. 

One might ask whether the ADM or harmonic coordinate system is more preferable, for the adequate physical 
treatment of the relativistic time delay and deflection of light rays in the field of gravitational waves emitted by a 
localized source. Our point of view is that the coordinate system should be chosen in such a way to be simultaneously 
both ADM and harmonic. The reason for this is that an observer who is at rest with respect to the ADM coordinate 
system does not feel the gravitational force caused by gravitational waves. This means that if the instantaneous 
gravitational field of the localized source may be neglected, the observer fixed with respect to the ADM system can be 
considered to be in free fall. Hence, no artificial force need be applied to the observer in order to keep him at rest at 
the fixed coordinate point. The motion of such an observer is described by the extremely simple equation x"^ = const 
and there is no need to account for kinematic effects associated with the observer's motion. All these advantages 
are lost in the "canonical" harmonic gauge. An observer fixed with respect to that coordinate system must be kept 
at a fixed coordinate point by some external force to prevent his motion under the influence of gravitational waves. 
The existence of such a force is unnatural from physical and astronomical points of view. On the other hand, the 
"canonical" harmonic gauge has the advantage of a much simpler integration of the equations of light propagation than 
the "canonical" ADM gauge. One can see that the "canonical" ADM metric coefhcients (B1)-(B2) contain functions 
which depend on time t only. As will be clear from the procedure of integration of equations of light propagation 
described in the next section such "instantaneous" functions of time do not permit explicit integration of each specific 
term (only after summing all terms is the explicit integration possible). Fortunately, the classes of ADM and harmonic 
coordinate systems overlap and, for this reason, we can substantially benefit by choosing a coordinate system that is 
simultaneously both ADM and harmonic. This allows us to proceed in the following way. First we integrate equations 
of light propagation in the harmonic gauge and then apply coordinate transformations (B13)-(B14) which transform 
the pure harmonic coordinate system to the ADM one without violation of the harmonic gauge conditions. This 
simplifies the treatment of observable effects drastically. 
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IV. METHOD OF INTEGRATION OF THE EQUATIONS OF MOTION 



A. Useful Relationships 

We introduce astronomical coordinates x = a;* = (x^, x^, a;^) corresponding to the plane of the sky of the observer 
and based on a triad of the unit vectors (Iq, Jo,Ko). The vector Kq points from the observer toward the deflector, 
and the vectors Iq and Jo lie in the plane of the sky, being orthogonal to vector Kq. The vector Iq is directed to 
the east, and Jq points towards the north celestial pole. The origin of the coordinate system is chosen to lie at the 
barycenter of the deflector which emits gravitational waves (see Figure 

Another reference frame based on a triad of the unit vectors (I, J, K) rotated with respect to vectors (Iq, Jq, Kq) is 
useful as well. The vector K points from the observer toward the source of light, and the vectors I and J lie in the 
plane of the sky, being orthogonal to vector K, which is different from the plane of the sky being orthogonal to vector 
Kq. This is because the "plane of the sky" is actually a sphere, and vectors K and Kq point in different directions. 
Mutual orientation of one triad with respect to another one is determined by the following equations 



Iq = Icosfi + J sinfi , (16) 

Jo = —I cos 9 sin 17 + J cos 9 cos 51 + K sin 9 , (17) 

Kq = Isinflsinfi — Jsin6'cosr2 + Kcos0 , (18) 

where rotational angles and 9 are constant. 



To integrate the equations of propagation of electromagnetic waves in curved space-time we must resort to an 
approximation method. In the Newtonian approximation, the unperturbed trajectory of the light ray is a straight 
line: 

x\t)^x'j^{t)^xl + k'{t~to), (19) 

where Iq is the instant of time of the photon emission from the point with spatial coordinates Xg, and fc' = k is a 
constant unit vector tangent to the unperturbed trajectory and directed from the point of emission to the point of 
observation of photon (the vector k « — K). In the Newtonian approximation, the coordinate speed of the photon 
= fc' and is considered to be constant. 

It is convenient to introduce a new independent parameter r along the photon's trajectory according to the rule 




r = k-x = t-to + k-xo, (20) 

where the dot symbol between two vectors denotes the Euclidean dot product of two vectors. The moment to of the 
signal's emission corresponds to the numerical value of the parameter tq = k • xq, and the moment t* of the closest 
approach of the unperturbed trajectory of the photon to the origin of the coordinate system corresponds to the value 
T = (note that tq < if the source of light is behind the localaized source of gravitational waves). Thus, we find 



T = t-t\ ro = to-t*. (21) 

The variable r is negative from the point of emission up to the point of the closest approach, and is positive 
otherwise. The differential identity dt = dr is valid and for this reason the integration along ray's path with respect 
to time t can be replaced by the integration with respect to parameter r. Using parameter r, the equation of the 
unperturbed trajectory of light ray can be represented as 

x'(t) =a;iv(T) =fcV + r, (22) 

and the distance, r, of the photon from the origin of coordinate system is given by 

r = rjv(r) = y/r^+d^, (23) 

where the length of the constant (for a chosen light ray) transverse vector ^ = k x (xq x k) = k x (x x k) is called the 
impact parameter of the unperturbed trajectrory of the light ray, d = |^|, and the symbol " x " between two vectors 
denotes the Euclidean cross product. It is worth emphasizing that the vector is directed from the origin of the 
coordinate system toward the point of the closest approach of the unperturbed path of light ray to that origin. The 
relations 
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r + T^ , ro+To = , (24) 

r - r r-o - To 

also hold, and they are useful for presenting the results of integration of the light ray equations in different form. In 
particular, if we assume the strong inequalities d <^r, and d <C ro to hold, then 

T = r- — + .., T-o = -ro + — + 25 

2r zro 

which clearly shows that at the moment of light reception t is positive and at that of light emission tq is negative. 

Let us consider a set of curves a;*(T) = fcV + ^' with different values of vectors fc* and The vector field fc% defined 
along the curve x^ir), describes the direction of a bundle of light rays along the curve, and introduces a natural 
"2+1" splitting of 3-dimensional space. The vector on the plane orthogonal to the bundle of light rays, is a point 
of intersection of any of those rays with that plane (see Figure |l|). This vector does not depend on r and can be 



defined, as in equation (22), by the relationship 

C = P}^' , (26) 

where 



(27) 



is the projection operator onto the plane orthogonal to the vector fc*. The operator has only two algebraically 
independent components and satisfies the relationship 



PlPf 



P- 



(28) 



Because of this property we can recast equation (pq) into the form 



r = p-e 



(29) 



which shows explicitly that the vector is constrained to lie in a 2-dimensional plane. Thus, we immediately have 
for the operation of partial differentiation in this plane 



^ _ pi _ pij _ p. . 



(30) 



It is worth noting that the projection operator can be used to raise and lower indices of any geometrical object lying 
in the plane orthogonal to vector fc' . 

In what follows, it is convenient to consider the spatial components of coordinates ^* as formally independent with 
subsequent projection onto the plane when doing differentiation with respect to Therefore we always use the 
operator of differentiation with respect to ^* in combination with the projection operator P^*. For example, before the 
projection we treat 



-5- = J' 



and for the same expression with subsequent projection 



pi _ pi 
3 g^q 3 



(31) 



(32) 



which agrees with equations (|2^) and 
function F(t, x) holds 



Moreover, the following rule of differentiation for an arbitrary smooth 
d , d 



x=xo+k(t — to) 



(33) 
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Equation ( |33| ) is a generalization of the corresponding formula introduced by Kopeikin ( [gij, equation (20)) for 
functions which do not depend explicitly on time t. It is worth noting that in the left-hand side of formula ( |33|) one 
has first to differentiate the function F(t,x.) with respect to time t and spatial coordinates and, then, to make the 
substitution x = xq + k(t — Iq). However, one makes corresponding substitutions in the right-hand side of the formula 
( ^ ) first and only afterwards takes derivatives. 

It is useful to stress again that because the coordinates f * lie in the plane orthogonal to the vector fc* only two of 
the three £,^,£,'^,£,^ are, in fact, independent. We also stress that the new variables and r are independent as well. 
For this reason, the integration of any function, which can be represented as a time derivative with respect to the 
parameter r, is always quite straightforward 

— F(r,Orfr = F(r,0 + C(0, (34) 

where C(^) is an arbitrary function of the constant impact parameter. Moreover, as the vector does not depend 
on time r, the partial derivatives with respect to ^* can be removed from within the time integrals when calculating 
them along the photon's trajectory, that is 

I^Fi-^i)dr^^jF{r,^)dr. (35) 

Because of these advantages the new independent coordinates r and are quite useful in calculations. The usefulness 
of the variables r and has been also recognized by Damour & Esposito-Farese |Q . 

The equations of motion of light rays (^ in terms of parameters ^ and r are simpler, and after accounting for a 
freedom in gauge transformations and implementation of relationship (^3|) assume the form p2| 

F(r) = ^^^h'o^■ - dr^r - Ik'drh'oT- " k'drh'ir + k'd,h'oT' + (36) 

i {^^ + k'dr) kn'ih''^^- - drr {w' - k' W°) , 



where the following notations are used: di = Pijd/dS,^ , dr = d/dr. Let us emphasize once again that the repre- 
sentation of equation (|3^) is valid in an arbitrary coordinate system and all metric coefficients are taken along the 
unperturbed trajectory of propagation of the light ray; that is, /iQ^(i,x) = hapir,^ -\- kr). We also remark that the 
right-hand side of equation ( |3^ ) contains only spatial partial derivatives with the same index "i" as does the left-hand 
side of the equation. This contrasts with equation (^ where the indices of spatial derivatives are mixed. Equation 
( ^ ) will be used in sections 5 and 6 for a general treatment of gravitational perturbations of the photon's trajectory 
and discussion of relativistic time delay and angle of light deflection. 

Another useful form of equation ( p6| ) may be obtained if one introduces the four- vector fc" — (1, /c*). Then we find 

i\T) = h^kPd.h'^^^- - dr (^k^Kr- - h'k^kPqf;''-^ ~ drr {w' - k' «^°) . (37) 

This form of the equation clearly shows that only the first term on the right-hand side contributes to the deflection 
of light, if observer and source of light are at infinity. Indeed, one integration of (^) with respect to time from —oo 
to -\-oo brings all first and second time derivatives to zero, due to the asymptotic flatness of the metric tensor. This 
makes a connection between the formalism of the present paper and that of Damour & Esposito-Farese (see also 
Appendix A). 



B. Calculation of Integrals from the Static Part of the Gravitational Field 

The static part of the gravitational field of the deflector contributes to perturbations of light's ray trajectory, defined 
by the following indefinite integrals [p3 



Mr, O^J'i-J = - In - r) , (38) 

B{t, €) = y Mt, ^dr = -T In ( v^cPT^ " ^) " Vd^+r^, (39) 
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where we have omitted constants of integration which are absorbed by re-definition of constants of integration of 
unperturbed hght trajectory (|l^). Integrals (^), ( ^9|) are formally divergent at the lower limit. However, this 
divergence is not dangerous for setting the second of the boundary conditions (|l]) because only derivatives of the 
integral (^) appear in the result of the first time integration of the equations of motion of light rays, eliminating the 
divergent part of the integral Esj. With this in mind, it is easy to prove that integrals (js^), ( p9| ) are in agreement 
with the boundary conditions 



C. Calculation of Integrals from Time Dependent Part of Gravitational Field 

One meets two ways of calculation of integrals in finding the path of propagation of light in the gravitational field of 



a localized source emitting gravitational waves. The first method relies upon the use of the Fourier transform (40) and 
allows one, at least in principle, to calculate all integrals explicitly if one knows the specific structure of the Fourier 
image of the quadrupole moment of the deflector [Q. The advantage of the second method is based on the fact 
that one deals with the metric depending on retarded time only. This allows one to make a special transformation 
of variables within the integral which excludes any dependence of the integrands on the impact parameter, and 
transfers it to the limits of the integrals. Thus, partial derivatives of the integrals can be calculated explicitly without 
assumptions about the structure of the quadrupole moment of the deflector. Of course, both methods give the same 
results. However, the second method is more general. 



1. First Method of Integration 



Let US assume the most general aperiodic form for the time variation of the deflector. In linear approximation the 
total mass and spin of the deflector are conserved quantities so that they do not depend on time at all, and we can 
consider them as contributing only to the static part of the gravitational field of the deflector |Q . The quadrupole 
moment is not static. It may be represented through a Fourier transform as 

/ + 00 
4(c^)e*"(*-'')du;, (40) 
-oo 

where Tij{oj) is the (complex) Fourier image of the quadrupole moment of the deflector which must be specified for 
any particular source of gravitational waves. Here, we need not know the specific structure of Tijiuj) as it will be 
shown later it is irrelevant for subsequent calculations. 
Taking time derivatives of the quadrupole moment yields 

/ + 00 
{iuj)%j{uj)e'^^'-'-Uu, (41) 



/ + 00 
{-uj'')%,{Lj)e'^'''-''^duj. (42) 
-oo 

Generally speaking, arbitrary aperiodic source of gravitational waves have an infinite spectrum. However, it is possible 
to choose that frequency band which gives the largest contribution to the spectrum. The mean frequency O of this band 
defines the size of far (wave) zone of the source, as being roughly equal to the wavelength of emitted gravitational 
waves A = 27rc/r2. For example, if the deflector of light rays is a binary system, then the strongest emission of 
gravitational waves takes place at twice the mean orbital frequency of the system. For making estimates we can use 
the following approximations for components of the quadrupole moment 



\r^^{Maec)j, ir^l ~ (X ec^) ^ , etc., (43) 

where a is a characteristic size of the source of gravitational waves and e is its oblateness, quantifying the deviation 
of the density distribution from spherical symmetry. 

When integrating the equations of light propagation using the metric with Fourier transform ( ^ ) for the quadrupole 
moment one meets the following integrals: 
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cos[a;(T — \/(P+T^)] 



dr. 



(44) 



l2iT,i,uj)= I . ^ dr. 



Vd2 + t2 



(45) 



In order to evaluate the integrals (p4|)-([45|) it is useful to change the time argument, r, to the argument y, by the 
transformation 



W2 



which yields 



2/2 - d 



1 d2 + j/2 

, V - . • : ai — „ 

2y ' 2 y ' 2 y2 



dr = ^ — dy. 



(46) 



(47) 



While the parameter r runs from — cx) to +oo, the new parameter y runs from — oo to 0; that is, y is always negative. 

After transforming time arguments, the integrals Ii and I2 are reduced to the cosine- and sine integrals respectively 
( formula 8.230): 



IiiT,^,Lu) = -Ciiuy), 



(48) 



I2{t,^,uj) = -Si(wy), 



(49) 



where constants of integration have been omitted. Secondary integration of integrals (p^)-(|49|) along the light trajec- 
tory is required as well. Using transformations (46)-(E^ we obtain 



Ji{t,^,uj)= I Ii{T,i,uj)dT^-TCi{ujy) + ^ujd^ 



Siiuy) 



cos{ujy) 

2y 



sin(^?/) 
2lo 



(50) 



J2(r,|,w)= / l2{T,i,u)dT^-T^i{ujy) + -ujd^ 



C\{ujy) 



sin([jy) 



2y 



-I- 



cos(wy) 
2uo ' 



(51) 



where constants of integration have again been omitted. 

Using the Fourier transform of the quadrupole moment (|4^) and formulae (^), (^), (^), ( ^9| ) one calculates the 
important integrals 



dt = (27r) 



/-I- 00 
4(t^)e^"** [h{T,^,u) + ih{T,tu:)]du, , 
-00 



(52) 



Q,(r,|) ^ r ^^2^Lj!ldt ^ (27r)-i/2 (^^^ oj%,{u)e^-'' [^h{T,£,,oj) + ih{T,£„uj)]doj , 



(53) 



/T /• + 00 

S,;,(T,Odi = (2^)-i/2 / 2-,,(^)e'-** [Ji(t, ^, c^) z J2(t, ^ c^)] d^ 
-00 J —00 



(54) 



/r /' + 00 

C,,(T,Odi - (27r)-i/2 / ^i^^.(^)e-** [J2(r, |, c^) - jJi(t, |, u;)] du; , 
-00 J — 00 



(55) 
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where t* is the moment of closest approach of the photon to the origin of coordinate system. In what follows, we need 
only partial derivatives with respect to the impact parameter of the integrals (|5^) - (|55|). These can be calculated 
rather easily. We have, for example, 

dMr,tio)^{yr)-^cosiLuy)e, a,/2(r, 1, ^) = (yr)-^ sin (u;y) f , (56) 

and so on. Thus, making use of the inverse Fourier transform we obtain 

^kB.,,iT,0^iyr^'I^At~r)e, (57) 
'5.i?.,(r,0=(l-^)^^^^^, (58) 
4c,,(T,|) = (^/r)-^4(^-r)e^ (59) 

9.a,(r,0=(l-n?^^^^, (60) 

Calculation of partial derivatives from integrals Dij{T,^) and Eij(T,^) may be done without difhculty in a similar 
fashion using equations (pO|)-(^. 

2. Second Method of Integration 



The second method also uses the substitutions (|4^), (^). The integrals (^ - (|53|) are brought into the form 

r?a(i^,<=_ri(£±o.c, ,61, 

J —OC ^ J — CO ^ 

r5i<i^<i*.-/"M;±ii.c. (62, 

One sees that the integrands of the integrals do not depend on the parameters d and r at all. They are present 
only in the upper limit of integration. Hence, the integrals (|6|), (|6^ ) are functions of the variable y only, that is 
-Sy (T, I) = Bij{y) and Cij(r, |) = Cij{y). Making use of the transformations (|^, (0), the integrals (||), (|5|) are 
reduced to the expressions 



r 1 /"^ fp R -(C) 

J — oo ^J— oo ^ J-oo S 



(63) 



C,,{T,^)dt^- a,(C)dC + -/ -^dC (64) 

Hence, the integrals Dij(r, ^), Eij^r,^) are also functions of the variable ?/ only. 

We stress once again that our formalism holds true for arbitrary dependence of the quadrupole moment of the 
localized source on time, and includes the case of sources which produce bursts of gravitational radiation, such as 
supernova explosions or coalescence of binary systems, as well as periodic systems. Indeed, suppose that the burst 
starts at the moment ti and terminates at the moment t2 ■ We assume for simplicity that before and after the burst the 
quadrupole moment of the source is identically zero. During the burst, the tensor function Tij{t) describes the time 
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dependence of the quadrupole moment. Then all formulae derived in this paper hold, if we describe the quadrupole 
moment of the source as a product of two Heaviside step functions with the tensor function !Fij{t). Thus, for any 
moment of time we write 



I,j{t) = H{t - h)H{t2 - t)T^j{t) 
where the Heaviside step function is defined as follows 

//(t-r) = {J 



\it>T, 
otherwise. 



(65) 



(66) 



Time derivatives of the quadrupole moment are calculated taking into account that H{t — T) = 5{t — T) is the Dirac 
delta-function, and 5{t — ti)Tij{ti) = 5{t — t-2)Tij{t-2) = 0. This yields 



(t) = H{t - h)H{t2 - t)T^J{t) 
and similar formulae for higher derivatives. 



(67) 



It is evident from the structure of integrals (k31^)-(64) that taking partial derivatives of any of the foregoing integrals 



is reduced to taking the partial derivative with respect to y. In particular, we obtain 



y 



-djy= {yr) Ipq{t - r)(^ , 



(68) 



which exactly coincides with the result ( |57| ) derived above using the inverse Fourier transform method. Second and 
third partial derivatives of the function Bij (r, ^) with respect to the impact parameter will be useful subsequently. 
They are calculated making use of formula (pq). This yields 



yr 



(69) 



and 



drjkBpq{T,i) = {yrY 



^2 ^2 



£,iPjk "^^kPj] '^CiS.j^k S^iCjCfc 



(70) 



(yr) 



^iPjk 'iS.kPij '^S.i^jS.k 3^i£.j£.k 



Jpq{t~r) 

ipq{t-r) 



ipq{t~r) 



3 -^Pl 



yr 

We note that the formulae of partial differentiation of Cij{T, ^) look the same as for Bij{T, ^) after taking into account 
the fact that the integral (^2|) depends on the first time derivative of the quadrupole moment. The derivatives of 
the functionals Eij{T,^) and Dij{T,^) can be obtained using relationships (|63|)-(|64) and derviatives of Bij{T,^) and 
Cij{T,^). For example. 



d,Dpq{T,^)=C' 



djkDpqir, ^) = ^dkBpq{T, + 



dtjkP>pq{T,£,) = - 



Bp,iT,^) 



-dC 



-OQ S 

P'' + ^ ) dkBpq{T, I) + P^^d,Bpq{T, i) + ^^d,kBpq{T, i) 



yr 



(71) 
(72) 
(73) 
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It is worth emphasizing that the third partial derivative of Dpq{T,^) does not include the integral Bpq{T,^) by itself, 
as might be expected, but only its first and second derivatives. Therefore, the third partial derivative of Dpq(r, ^) 
does not depend on the past history of propagation of the light ray (see formulae ( |6^ ) and (|69|)). 

Now, after making these remarks, we are ready to discuss the relativistic perturbations of the photon's trajectory 
in the radiative gravitational field of a localized source deflecting light rays. 



V. PERTURBATIONS OF PHOTON'S TRAJECTORY 



We first note that in terms of the new variables r and the components of the "canonical" metric tensor (|T^)-(|TJ) 
taken at an arbitrary point on the light ray can be re- written as follows p8| : 



'00 



r,0 = 



r^can. , 



2M 

r 



ki kj df 



Jij (t-r 



{t - r) 



2 [di+kidr 



i%j{t - r) 



%j{t - r) 



2k 



(74) 



(75) 



TA)=5.,hlZ{rA) + -%,{t-r), 



(76) 



where in the right-hand side of all formulae it is implicitly assumed that variables i, are replaced by r and (J' , and 
di = Pldjd^^ dr = d/dr. In addition, note that the dot over the quadrupole moment takes the usual meaning of 
differentiation with respect to time, which must be completed first, before substitution of t and for r and and 
before taking any other derivative. 

The metric tensor (jfl) - (||) is used in the equations of motion of light rays (^6|) which are reduced with the help 
of formula (pSh to the expression: 



2M (di-kidr) -2SP 



kqCipqdfT kj CjpqO-iq 



(77) 



(^dipq kidpq-T -\- '2kpOiqj- ~h kpkqdifr '^kikpQqf^ kikpkqd-j-j-j^ 



^ki kpdqf 



^ip^QT ^ipkqdi-r H~ kikpkqdff 



lpq{t-r) 



where and w° are functions given by relationships (|^)-(10). Remarkably, no terms depending on the second 
time derivatives of the quadrupole moment appear in the equations of motion of light rays ( [TT] ) , because of mutual 
cancellation. This fact explicitly demonstrates that gravitational waves emitted by localized sources are much more 
elusive from detection by angular deflection than other authors suggest. It is worth noting that the disappearance of 
terms with second derivatives from the quadrupole moment is a local phenomena and is not a result of integration 
of equation ( |77| ) . This is a characteristic feature of General Relativity. Alternative theories of gravity do not possess 
such a local cancellation of gravitational wave terms. This cancellation may be used for conducting new tests of 
General Relativity in the weak, radiative grav itational- field limit. 

Let us simplify the equations of motion (|77|) in order to avoid writing down cumbersome expressions. We introduce 
two functions and (f^^ which generate in (|77|) the time derivatives of second and higher orders. These functions are 
defined: 



(^0 = - 2A:pV, 



= 2SP kg Er 



-\-kpkq 



r\jr> rxjQ V i 



^kn 



itq{t - r) 



(78) 
(79) 
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where the differential operator Vi = d/dx^ must be apphed before the substitution of the unperturbed trajectory of 
hght rays. It can be easily confirmed by straightforward use of formula (|3|) that the expressions (|78|)-(|79|) generate 
terms with second and third derivatives with respect to r in (^^. The equations for the path of the light ray now 
assume the form: 



(80) 



drr [w' +ip'-k' + 



We note that the terms ip'^ and ip^ are gauge-dependent and can be, in principle, eliminated from the equations of 
motion (|77| ) by choosing appropriate gauge functions and w*. However, such a procedure will introduce a reference 
system with a coordinate grid very sensitive to the direction to a specific source of light rays; that is, to the vector A:*. 
The coordinate system obtained in this way will be of trifling practical usage. For this reason we do not recommend 
the elimination of functions ip^ and (p^ from (77) and give preference to the ADM- harmonic coordinate system, which 
admits a much simpler and unique treatment of observable effects. Thus, we leave the functions ip'^ and in the 
equations of motion of light rays, where gauge functions and are defined by formulae ( B13| )-(B14). 
Proceeding further in this way and integrating equations (^7|) one obtains 



(81) 



x\T)=xUr) + E\T)-E\To) 



(82) 



where the unperturbed trajectory of light ray x]^{t) is determined by the expression (|22|). The relativistic perturba- 
tions to the trajectory are: 



E\t) = [2Md, + 2SPkjejpqd,qj A{t,^) + d,pqBpq{T,^) 

(2Mh + 2SPe,pqdq) 



dr [w'+^'-k' {w°+ip°)] 



Ipqit-r) 



-2P 8 



(83) 



ijqit-r) 



S'(r) 



2Md, + 2SPk,ejpqd^q) B{t,^) - (2Mh - S^^'e.p,^) A(t,|) 

dipqDpq{T,C) - {kidpq - 2kpdi^ Bpq{T,i) - 2PijdqCjq{T, $,) - 



(84) 



We emphasize that before differentiation with respect to time r or impact parameter 
quadrupole moment with respect to time t and make the substitutions: i i— > r, r i-^- Vd^ 



one has to diff erentiat e the 
d^+T^. We 
|84|) are A(r, ^) 



also wish to underline that the only integrals which need be calculated explicitly in expressions 
and B{t,^). All other integrals are acted upon by partial derivatives, which reduce them to ordinary functions as 
explained in the previous section. This remarkable fact allows considerable simplification of the calculations. This 
simplification results from the fact that the integrands can be formed from retarded potentials independent of impact 
parameter, after using the transformation of variables (^). This would be impossible if the metric tensor were not a 
function of retarded time t—r. Thus, retardation simplifies the calculations in the case of time-dependent gravitational 
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fields. In the case of a static or stationary gravitational field, the calculation of propagation of light can be done using 
the same technique since one can always consider a constant multipole also as a (constant) function of retarded time. 
For this reason, more involved calculations of light propagation (e.g. see M and pa]) can be simplified as well. 



Th e functions w"^ and w^, which describe freedom in choosing coordinate systems, are taken from formulae ( B13 ) 
(B14) of Appendix B. Consequently, the integrals of equations of light propagation ( [77| ) expressed in the ADM 



harmonic coordinate gauge possess a simple interpretation of observable effects, as discussed in the following section. 

It is convenient to obtain an expression for unit vector fc' written in terms of spatial coordinates of the points of 
emission, Xq, and observation, x, of the light ray. From formula (|8^) one has 

|x- Xo| 

or more explicitly 

k'^~K'~/3\r,^)+P\To,i) , (86) 
P'ir, = /3m (r, + Phir, I) + /3^(r, I) , (87) 

where the relativistic corrections /3*(t, ^) to the vector are defined as follows: 

a,(r.O.?^fi^. (88) 

|X Xq I 

^ ' |x-xo| 

_ di^Mr, + 2kpd,,Bp,{T, ~ 2P^W,C,,{t, |) - Pj p (T, ^) + ^^ (r, |)] 

[r,i)- ^^3^^^ . (90) 

The relativistic corrections /3'(to,^) are obtained by replacing the parameter r in the numerators of expressions 
(^)-(|90|) by tq. One notes that in equation ( ^6| ) the unit Euclidean vector 

if' = -f^ (91) 
|x-xo| 

defines the direction from the observer towards the source of light and may be interpreted as a direction in asymp- 



totically flat space-time Relationship (85) allows us to apply the results of integration of equation of light 

propagation to the boundary value problem as well. The boundary value problem is formulated in terms of initial xq 
and final x positions of the photon 

x(t) = X , x(tn) = xo , (92) 

whilst the initial-boundary value problem (|l|) is formulated by means of assignment of the initial position xq and 
velocity of photon at past null infinity. The relativistic correction to the vector K"^ contains in its denominator the 
large numerical value of the distance between observer and source of light. However, the difference 'BP (r, ^) — (tq, ^) 
in the numerator of (|8^) may be of the same order as |x — Xo| itself. For this reason the relativistic correction in 
question must be taken into account, in general, for calculation of light deflection in the cases of finite distances of 
observer or source of light from the localized source of gravitational waves. Only in the case where observer and 
source of light reside at large distances on opposite sides of the source of gravitational waves, as was assumed in the 
paper by Damour & Esposito-Farese (1998), can the relativistic correction /3* beneglected. 

VI. BASIC OBSERVABLE RELATIVISTIC EFFECTS 
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A. Time Delay 



The gravitational time delay is derived from equations (^2|), (^). In order to obtain the expression for the time 
delay we multiply equation ( |82| ) by itself and then find the difference t — to hy taking the square root and using an 
expansion with respect to small relativistic parameters. This yields: 



i - io = |x - xol - k • H(t) + k • H(ro) 



(93) 



or 



t - to = |x - xol + AM(t, to) + As{t, to) + Aqit, to), 



(94) 



where |x — xo| is the usual Euclidean distance | |50|] between the points of emission, xq, and reception, x, of the 
photon. Am is the classical Shapiro delay produced by the (constant) spherically symmetric part of the gravitational 
field of the deflector, is the Lense-Thirring or Kerr delay due to the (constant) spin of the localized source of 
gravitational waves, and Aq describes an additional delay caused by the time dependent quadrupole moment of the 
source. Specifically we obtain: 



Am = 2M In 



ro + To 



(95) 



As = -2eyfcF5'=a, In 



ro + To 



(96) 



where 



Aq - % [B,, (r, I) - B,, (ro, ^)] + Sq{t, |) - 6q{to, ^) 



Sq{t, $,) = (w' + ifi') - - = 
^-'%,it~r) 



(97) 

(98) 



(-%,(< -r) 



and functions ^^Xpq{t — r) and ^ '^^'Ipq{t — r) are defined by formula (B5) of Appendix B. The expression for the 



second derivative of function Bij{T, ^) has been given in equation (|69|). The other derivatives appearing in Aq are as 
follows 



dr 



I,j {t - r) 



y ljjjt - r) _ T J,j{t - r) 



(99) 



^-^%q{t-r) 



Ipqit - r) + 3' 



Ipq{t-r) ^ (-i)Tp,(t-r) 



(100) 
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VpV, 



,(-i)Xp,(i-r) , (-2)j^,(t-r) 



P Q 



(101) 



dr{ VpV, 



Ipq - r) + 3 



-3^ 
r 



ipq{t - r) + 3' 



(102) 



The relationship (^4|) for the time delay has been derived with respect to coordinate time t. In order to convert 
this relationship to observable proper time, we assume for simplicity that the observer is in a state of free fall and 
that his velocity is negligibly small at the point of observation, with spatial coordinate x. If the observer's velocity 
is not small an additional Lorentz transformation of time must be applied. Transformation from the ADM-harmonic 
coordinate time t to proper time T is made with the help of the formula (e.g. see ||51|) 



(103) 



Implementation of formula ( |B7D for /iqo and subsequent integration of (|103D with respect to time yields 



T 



U) 



(104) 



where t\ is the initial epoch of observation and all velocity-dependent terms are assumed small, as argued above, and 
are therefore omitted. We also stress that under usual circumstances the distance r is so large that the difference 
between the observer's proper time and coordinate time can be neglected. Thus, we are allowed to treat coordinate 
time t as proper time. 

We note that the time delay in the propagation of light depends not only on instantaneous functions of retarded 
time but also on the integrals of time ''~^-'Ipg(i — r) and '■^^^Xpq{t — r). These integrals describe the whole past history 
of the source of gravitational waves up to the moment of observation. Under usual circumstances, the influence of 
such integrals on the time delay is expected to be small. However, this question deserves more detailed discussion and 
will be studied in more detail elsewhere. For example, these terms may be revealed in observations as the "kinematic 
resonance effect" predicted by Braginsky & Grishchuk |]l^. These terms may be also important for detection of the 
solar g-mode tidal oscillations by the LISA gravitational- wave antenna in space p^ . 

B. Deflection of Light 

The coordinate direction to the source of light measured at the point of observation x is defined by the four- vector 
— (l,p*) where ~ — i*, or 



P 



(105) 



and the minus sign directs the vector from observer to the source of light. However, the coordinate direction p^ is 
not a directly observable quantity. A real observable vector towards the source of light, = (1, s*), is defined with 
respect to the local inertial frame of the observer. In this frame = —dX'^/dT, where T is the observer's proper time 
and X"^ are spatial coordinates of the local inertial frame. We shall assume for simplicity that observer is at rest |5^] 
with respect to the (global) ADM-harmonic coordinate system (<, x*). Then the infinitesimal transformation from the 
global ADM-harmonic coordinates (i,x*) to the local coordinates (T, X*) is given by the formulas 



dT = Kldt + A° dx^ 



dX' = Aidt + a; dx^ 



(106) 
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where the matrix of transformation is defined by the requirements of orthonormality 



(107) 



In particular, the orthonormality condition (107) assumes that spatial angles and lengths at the point of observations 
are measured with the Euclidean metric 5ij. Because the vector s" is isotropic, we conclude that the Euclidean length 
|s| of the vector is equal to 1. Indeed, one has 



-1 + = 



(108) 



Hence, |s| = 1. 

In the linear approximation with respect to G, the matrix of the transformation is as follows |31 



Ao = 1 - 2'^oo(i,x) , 

A? --/io.(i,x) , 
A^ = 0, 



A'' = 

3 



1 + 2^oo(i,x) 



5^, + ^hl^{t.^) 



(109) 



Using the transformation (106) we obtain the relationship between the observable vector s* and the coordinate direction 



In the linear approximation this takes the form 



(110) 



s' = {l + hoo~hoy)p' + hj;^p^ . 
Remembering that vector |s| = 1, we find the Euclidean norm of the vector from the relationship 



(111) 



IpI = 1 - /loo + hojP" - ^hj^p'pi 



which brings equation (111) to the form 



s* = TO' + ip*-'m«/iJ/^(t,x) , 



(112) 



(113) 



where the Euclidean unit vector = pY|p|. 

Let us now denote by a* the dimensionless vector describing the total angle of deflection of the light ray measured 
at the point of observation, and calculated with respect to vector fc' given at past null infinity. It is defined according 
to the relationship 132] 



a^(r,0 = fc^[k-H(r,0]-SXr,C) 

or 

a^(r,0--/^*S^(r,0- 
As a consequence of the definitions ( |105| ) and ( |115| ) we conclude that 



(114) 
(115) 
(116) 
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Accounting for expressions (113), (116), and ( ^6| ) we obtain for the observed direction to the source of hght 

s^(t, ^) = + a^T, + I3\r, ^) - /3\to, + ^\t, , (117) 
where relativistic corrections (3^ are defined by equations (|8l)-(p0|) and the perturbation 



7 



{TA) = -\p''k'^h^^{t,^) 



(118) 



If two sources of hght (quasars) are observed along the directions s\ and s\ the measured angle ip between them in 
the local inertial frame is: 



cos ^/j = Si • S2 , 



(119) 



where the dot denotes the usual Euclidean scalar product. It is worth emphasizing that the observed direction to the 
source of light includes the relativistic deflection of the light ray. This depends not only on quantities at the point 
of ob servation but also on /3*(ro,^), at the point of emission of light. This remark reveals that according to relation 
(119) a single gravitational wave signal may cause different angular displacements for different sources of light located 
at different distances from the source of gravitational waves. 

Without going into further details of the observational procedure we give an explicit expression for the angle a*. 
We have 



a" a\i (r, I ) + a's (r, ^ ) + (r, ^ ) , 



(120) 



where 



ali{T,^)^-2Md^AiT,^) 



"q('^'^) = ~ dtpqBpq{T,$,) - Pij [2kpdjq + kpkqdjr + 2kpSjqdrr + 25jqdp 



(-2)j,,(t-r) 



(121) 

(122) 
(123) 



The expression for the third spatial derivative of function Bpq (r, ^) has been given in equation ( [70| ) . The other relevant 
derivatives are: 



(124) 



dr 



y Ijj {t - r) _ T Ijj (t - r) 



(125) 



Ipq{t-r) 



iq 



ipq{t - r) Ipq{t ~ r) 



(126) 



^i^q 



ipq{t-r) 3ipq{t-r) 3Ipq{t-r) 
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ipqit ~ r) ipqit - r) 



(127) 



dr7 



(128) 



Straightforward but tedious calculation of the last term in equation (123) yields 



9^ VpV,-"'^^^^(* 



1pq{t - r) + 3 



(-i)lp,(t-r) , (-2)lp,(t-r) 



P 



Ipqit - r) +6 



Ipg{t-r) , ^J-^%q{t-r) , ^(-^%g{t-r) 



15- 



15- 



fri 



(129) 



and 



d,r\ VpV, 



i-'%q(t-r) \ 



= 2 



-6^ 

r 



Ipqit - r) + 3 
Ipqit~r) + 4: 



^-'%qit-r) ^ , J-'%qit-ry 
^-'%qit-r) ^ ^ ^-'%qit-r) 



kpp, 



iq 



P tP 
^ iq-^ TV 



-2^ 
r 



r 



(-i)2p^(t-r) 



r 

kPP 



iq 



-2 



ipg(i - r) + 6 



^pqjt - r) ^ ^^ (-^%qit-r) ^ ^^ ('^%^{t-r) ^ Ck^xj, 



+4- 



15 Jpgit-r) 45 (-i)Tp,(t-r) 45 (-^)Xp,(t - r) ' 
2 r 2 r2 2 



.P ™9 



ipqit - r) + 6 



(130) 



We note that the angular displacement in astrometric positions of sources of light in the sky depends not only 
on quantities that are instantaneous functions of retarded time, but also on integrals over time '^"^^Ipqit — r) and 
^'"^^Xpqit — r), which describe the whole past history of the source of gravitational waves up to the moment of 
observation. Under usual circumstances the influence of such integrals on the deflection of light is expected to be 
small. However, this question deserves more detailed discussion and will be discussed elsewhere. 



VII. DISCUSSION 



It is remarkable that amo ng aU the integrals required for calculation of the trajectory of the light ray, only i?y (r, ^) 
enters the expressions (|9^), (123) for time delay and deflection angle. Furthermore, it is remarkable that we need not 
know this integral explicitly, but only its second and third derivatives with respect to impact parameter. These are 
given in equations (69) and (^3)- With the knowledge of these derivatives, and derivatives of other functions given in 
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the previous section, we have complete information about the functional structure of the relativistic time delay and 
the angle of light deflection produced by any localized gravitating system possessing a time-dependent quadrupole 
moment Tij (t) . 

This structure indicates that the explicit time dependence of the quadrupole moment completely determines the 
results of astrometric and timing observations. We shall not consider this problem in the present paper, leaving it for 
future exploration. 

Our concern in this section is the simplification of the general formalism developed in the foregoing text. In order 
to do this we consider three limiting cases: 

1. The impact parameter d is much smaller than the distance from the localized source of gravitational waves to 
both the observer, r, and to the source of light, vq. The source of light is behind the source of gravitational 
waves (see Figure ||); 

2. The impact parameter d is much smaller than the distance from the localized source of gravitational waves to the 
observer, r, and to the source of light, rg. The source of light is on the same side of the source of gravitational 
waves as the observer (see Figure ||) ; 

3. The distance R from the source of light rays to the observer is much smaller than distances from the observer 
or from the source of light to the localized source of gravitational waves. The impact parameter d may be 
comparable with the distance from the deflector to observer or the source of light (see Figure ^) . 

We will conventionally refer to the cases 1 and 2 as those of small impact parameter, with numerical values of tq < 
and Tq > respectively. Case 3 is that of large impact parameter, and its small numerical values are covered by the 
formalism as well, as will be clear in section 7.3 below. 

A. Case 1. Small Impact Parameter (ro < 0) 

1. Asymptotic expansions of independent variables 

We shall assume in this section that the condition d <^ min[r, tq] holds. Let L = min[r, ro] and recall that 
T = y/r^~~(P and tq = — \/r^~~(P < (see Figure ||) . This yields 

7^ d'^ d'^ , , 

y^^r^-d^-r^-— - — + ..., 131 
2r 8r^ 



and 



yo = -\lrl -d^-ro = -2ro + ^ + ^ + - . (132) 



where dots denote terms of higher order, r is the constant distance from the deflector to observer, and ro is the 
constant distance from the deflector to the point of emission of light. Using these expansions we find: 

t = r+r - — + ..., io=r-ro + — + ... . 133 

2r 2ro 

These can be used for Taylor expansion of functions about the time t* , the moment of the closest approach of light 
ray to the deflector. Specifically, if we assume convergence of this Taylor series we find: 

J.,(t - r) = J.,(i*) - ^i,(r) + ... , (134) 



2ro' 



(to - ro) = J,, (r - 2ro) + —X,, {t* - 2ro) + ... , (135) 



23 



where dots again denote terms of higher order. Convergence of the time series given above requires: 



« 1 , 



and 



c ro 



«1, 



(136) 



where uj is the highest frequency of gravitational waves emitted by t he d eflector. If the source of light rays and observer 
are at infinite distances from the deflector then the requirements ( 136 ) are satisfied automatically, irrespective of the 
structure of the Fourier spectrum of the quadrupole moment of the deflector. In practical situations such an assumption 
may not be always satisfied. For this reason, it will be more natural to avoid the Taylor expansions of the quadrupole 
moment with respect to retarded time. It is also worth noting that in the case of small impact parameter we have 



[yry 



and 



1 i 



2rl 



8r4 



The foregoing expansions then yield 



(137) 



(138) 



djBpgir, = -2d, In + 1^ j I,fe(i - r) 



2rl 



1jk{to - fo) + 



djkBpqir, I) = -2djk lndJpq{t - r) + -ujUkipqit - r) 



djkBpq{To,$,) — ——^Pjklpq{ta — ro) 



di-jkBpq[T,Ci 



dijkBpq{To,^) =0 -3 , 



^pq{t - r) + —Ipqit - r) 



1 



dijk Ind 



(139) 
(140) 
(141) 
(142) 
(143) 
(144) 



and 



dijkDpq{T,i) = -2rlpq{t - r)di-jk Ind 



dijkDpq{To,i) = O ( — 



1pq{t - r) + 



(145) 
(146) 



In addition we have 

5Q{T,(,) = -kn'>±pq{t-r) + ..., (147) 

5q(to,0 = o(^^) . (148) 

We note that the leading terms of the expansions decay much faster (at least as I/t-q) at the point of emission of light 
than those at the point of observation. This indicates that the main contribution to the effects of time delay and 
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deflection of light arise along the path of the light ray from the localized source of gravitational waves to the observer. 
We discuss this question in more detail in the following section. 

The asymptotic expansions of integrals (|3|) - (|9|) describing propagation of light rays in the static part of gravi- 
tational field of the deflector are: 



A{t,^) = -21nd + ln2r 



4^2 



(149) 



4r2 



(150) 



Bir,^) = -r-2rlnd + rln2r- 



2r 



2 \2r 



(151) 



/I 

B{tq , ^ ) = -ro + ro In 2ro - — ( - + In 2ro 



(152) 



These expansions are used for calculation of asymptotic expressions for time delay and the angle of deflection of light 
rays. 

2. Asymptotic expressions for time delay and the angle of light deflection 
The static part of time delay and deflection angle are: 



Am = —iM liid+2M ln(4rro) 



As = -4e,,pP5Pa 



Ind — — ln(4rro) 



In d — — ln(4rro) 



In d — — ln(4rro) 



(153) 
(154) 
(155) 

(156) 
(157) 
(158) 



where we have neglected the angle 7*(t, ^) because it is small (recall that 7'(r, ^) ~ P^^ k'^hJJ'). 

Asymptotic expressions for the time delay and angle of deflection caused by the quadrupole moment are: 



A, 



'22ij{t — r)dij Ind H — {2ninj + kikj)2ij{t ~ r) + ... , 



(159) 



and 



25 



^Mt - r) + -Ijk{t - r) 



dijk In d 



k'ljk [t " r)dik In d + -C^jk [t ~ r)djk In d 



(160) 
(161) 



where — ^ /d is the unit vector directed along the impact parameter, R — |x — xo|, and dots denote terms of 
higher order [Q. The angle /3'(to, ^) at the point of emission of light is negligibly small and, for this reason, its exact 
expression has been not shown. 

Our calculations show that the time dependent part of the time delay and light deflection by the quadrupole moment 
of a localized source of gravitational field fall off in the first approximation as the inverse square and inverse cube 
of the impact parameter d respectively. For this reason there is no magnification of the gravitational wave signal in 
astrometric or pulsar timing observations as some authors have suggested pj|] - [^l|. In particular, terms proportional 
to l/d, or even to 1/d^ , appear only in terms of high order in the expansion (160) and are always multiplied by the 
factor 1/r to some power. 

The first term of formula (159) was first derived by Sazhin 0] for the special case of a binary system with a specific 
orientation of its orbital plane. Our derivation of formula (160) improves and gives independent confirmation of 
the result established previously by Damour & Esposito-Farese p3[ using another mathematical technique based on 
application of Fourier transform and pure harmonic coordinates. For completeness we have repeated the calculations 
of Damour & Esposito-Farese (23) for the effect of deflection of light rays by localized sources of gravitational waves 
in ADM rather than harmonic coordinates (see Appendix A) . The result coincides completely with that of Damour 
& Esposito-Farese |23 and clearly demonstrates the gauge invariance of the result. However, our technique is more 
general and powerful. Our formalism is valid for any relative position of observer, source of light, and source of 
gravitational waves, and with finite or infinite separations. The method developed by Damour & Esposito-Farese 
p3[ is valid only for infinite separations and for small values of impact parameter. In particular, we note that while 
Damour & Esposito-Farese find that the deflection depends on the time t* of the closest approach of light to 
the deflector, our calculation shows that it depends on the retarded time t — r. This difference is insigniflcant for 
extremely large separation of the light source and observer from the deflector, and small impact parameter, but it 
can be important in the cases of flnite distances or large impact parameter. 

It is important to realize that in the case of a small impact parameter, the basic time-dependent contribution to 
the bending of light and time delay by the gravitational field of a localized source of gravitational waves comes from 
the static part of the near-zone gravitational field of the source taken at the retarded time (cf. formulae (50)- (53) 
from ||2^ ) . In this respect it is worth emphasizing that the formula for the bending of light given in paper [^3| as well 
as in Appendix A is valid under two assumptions: 1) the impact parameter d is small compared with the distance to 
the observer r, 2) the velocity of matter inside the source of gravitational radiation is much smaller than the speed of 
light (the slow- motion approximation). 

The first assumption is rather trivial, since the impact parameter d is the only finite distance when the source of 
light and observer are at infinity. The second assumption appears because paper uses the Taylor expansion of 
the Fourier image of the tensor of energy-momentum of matter with respect to wave vector k (see equations (3.3) 
and (3.4) of paper p^). This expansion is mathematically equivalent to the use of a slow-motion approximation 
psf which, in particular, restricts the nature of the source of gravitational waves so that its Fourier spectrum is not 
allowed to include too high frequencies. 



In contrast, the general formalism given in the present paper produces results (159) and ( |160| ) applicable to arbitrary 
sources of gravitational waves, including gravitational radiation bursts with internal velocity of matter comparable to 
the speed of light [Q. Moreover, we do not assume positions of observer and the source of light to be at infinity. 

If we introduce the notion of the transverse-traceless (TT) and longitudinal (L) tensors pj], |Q with respect to 
the direction of propagation of light rays 



(^Sipkjkq -|- ^jpkikq^ ^p, 



(162) 



It 



k^kj {kpk(^p. 



(163) 



the expressions (|159D-(|160|) arc reduced to the form 



Aq = -2lf^{t - r) d,j lnd+ -n^nP ij^ {t - r) + 



(164) 
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dijk In d 



R 



F Jj^, (t - r)d,u In d + ^ - r)a,fe In d 



(165) 
(166) 



which reveals expHcitly that only the transverse-traceless part of the quadrupole moment of the localized source of 
gravitational waves contributes to the leading terms. However, terms of higher order can depend on the longitudinal 
part of the quadrupole moment as well. 



It is interesting to see that if we apply the expansions (134)-( |135| ), use the approximation of a gravitational lens, 
and omit all terms depending on time derivatives of the quadrupole moment, the expressions for the time delay and 
the angle of light deflection can be reduced to the formulae 



t - to = |x - xol - ^ + 27W ln(4rro) , 



(167) 



where '0 is the gravitational lens potential [p8[ having the form 



M 



in A 



pq 



\nd. 



(168) 



Scrutiny of the multipole structure of tp in cosmological gravitational lenses may reveal the presence of dark matter 
in the lens and identify the position of its center of mass, velocity and density distribution. 



Expression (168) includes explicit dependence on mass, spin, and quadrupole moment of the deflector and generalizes 
that given by Damour & Esposito-Farese |Q by accounting for the spin multipole. A similar result for the gravitational 
lens potential was obtained independently by Kopeikin ||2^ in the case of a stationary gravitational field for the 
deflector. The fact that the deflection angle can be represented as a gradient of the gravitational lens potential 
-0 explicitly indicates that the, so-called, frame-dragging effect in gravitational lenses can give a noticeable 
contribution to the deflection angle. Frame-dragging also produces a small displacement of the image of the background 
source from the plane formed by the two vectors directed from the observer toward the image of the light source and 
toward the gravitational lens. This torsional displacement of the image is pr oduc ed only by the component of spin 
of the deflector directed along the light ray (see the second term in equation (158). The overall effect of the torsion 
is of order d/r smaller than the main terms in the expression ( |168| ). These remarks dispel a seemingly common 
opinion that rotation of the deflector prevents representation of the deflection angle as a gradient of a gravitational 
lens potential. Similar conclusions can be derived from pi] and p3]. Ibahez fc Martin [|60| and Ibahez pll] give a 



formula for effects of frame-dragging equivalent to the spin-dependent term in (168), although they do not calculate 
all necessary integrals or estimate residual terms. 

Taking into account formula (117) and expressions for a% and 7 we obtain the vector equation for a gravitational 
lens 



R 



(169) 



where a' is given by relationships (167), (168) and we have taken into account that in the case under consideration 
i? ~ r -I- tq. One recognizes that when distances are finite the deflection angle with respect to vector is not simply 
but the product of tq/R and a*. In the limit when fc% which is equaivalent to — > 0, or r = const., tq — > 00 

the observed angle of deflection approaches the total angle of deflection a*, as it must in this limiting case. 



B. Case 2. Small Impact Parameter {tq > 0) 

1. Asymptotic expansions of independent variables 



We shall again assume in this section that the condition d ^ min[r, tq] holds and that r = \/r^ — d'^ and tq 
yjr\ - d^ > (see Figure |). This yields 
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(170) 



and 



yo 



2ro 8rg 



(171) 



where dots denote terms of higher order, r is the constant distance from the deflector to observer, and is the 
constant distance from the deflector to the point of emission of hght. Using these expansions we obtain the following 
decompositions 



t = t* + r- 



2r 



to = t* 



2ro 



(172) 



These can be used for Taylor expansion of functions that depend on retarded time about the time t* . In this case t* 
is the moment of closest approach of the light ray trajectory extrapolated backward to the deflector (see Figure ||). 
If we assume convergence of this Taylor series we find: 



(173) 



I,,(to-ro)=X,,(r)-^4(r) + ... 
where dots again denote terms of higher order. We also have 

, ,_i 2 1 

and 



(174) 



(175) 



iVoro) ^ 



(F 2rl Sri 



The foregoing expansions yield 



(176) 



djkBjk{T,$,) = -2djk Indljkit - r) + ^fijHkljkit - r) + 
djkBjkiro,^) 
dijkBjk{T,$,) = -2 
dijkBjk{To,C) — —2 



2dj In d + — ^ ) Xjkih - ro) + 



2djk \ndIjk{to - ro) + —njnkIjk{to - ^o) + 



d^ ■ 

Ijk{t - r) + —Ijk{t - r) 



dijk Ind- 



dijk In d - 



(177) 
(178) 
(179) 
(180) 
(181) 
(182) 
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and 



dijkDpq{T,^) = -2rlpg(t - r)dijk Ind 



ipq{t-r) 



dijkDpqijQ, ^) = -2roIpq{to - ro)i9yfe Inrf 1 ^pq(to - ^o) + 



(183) 
(184) 



In addition we have 



6Q{T,^)^-kPk'^Ipq{t~r) 



Sq{to,0 = -k^'k^Ipgito - ro) + 



(185) 
(186) 



We note that the leading terms of the expansions now have the same dependence on the distance of the point of 
emission of hght and of the point of observation from the source of gravitational waves. If the source of light is closer 
to the source of gravitational waves than the observer, it makes the largest contribution to the effects of time delay 
and deflection of light. 

The asymptotic expansions of integrals (|8|) - ( ^9|) describing propagation of light rays in the static part of the 
gravitational field of the deflector are: 



^(T,|) = -21nd + ln(2r)- — 



(187) 



^(ro , = -2 In d + ln(2ro) - 



(188) 



B{t, ^) = -r - 2r In d + r ln(2r) 



j2 r 



2r 
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(189) 



(190) 



These expansions are used for calculation of asymptotic expressions for time delay and the angle of deflection of light 
rays. 

2. Asymptotic expressions for time delay and the angle of light deflection 
The static part of time delay and deflection angle are given by: 



2M 



In 



4 



1 1 



(191) 



As = e.,pyS^C ( 4 - ;i 



(192) 
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Expressions for a\,j, a^, and Oq will be the same as in equations ( |155D , ( |l5q ), and ( |16C| ) because they are taken at the 
point of observation only. Expressions for j3\j^ Pg, and /3q are given at the point of observation by equations ( 157 ), 
(158), and (161). Expressions for (3\^, /3g, and /3q at the point of emission of light are given by the same equations 
(157), (158), and (161) after substituting ro for r. The relativistic perturbation 7* is calculated in equation (118). 
The asymptotic expression for the time delay caused by the quadrupolc moment is: 



Aq = -2 [li-j (t - r) - (to - ro)] 9^ In d 



(193) 



(2ninj 



ki kj ) 



ro 



(194) 



One might think that the effect of retardation is again inversely proportional to the square of impact parameter d. 
However, this is actually true only for sources of gravitational waves with rapidly varying quadrupolc moment. If 
motion of matter inside the localized source of gravitational waves is slow, then conditions (136) apply. In this case, 
the real amplitude of the effect becomes extremely small, being inversely proportional to l/r^ and l/rg. 

Th e asymptotic expression for the observed direction s* to the source of light is derived from the basic formula 
(ITTtI) and is: 



s^^K'--A [x,k{t - r) - I,k{to - ro)] d^,k In d 
H 



(195) 



[X,7c {t - r) - Ijk {to - ro)] dik In d 

ijk{t - r) - ijk{to - ro) djk \nd + f{T, |) + 



where we have accounted for the approximate equality i? ~ r — rg valid in the case of tq > 0. One can see that 
deflection angle is small in the expression given. Moreover, if we again assume that motion of the matter is slow, then 
the observed deflection is even smaller and is inversely proportional to l/(ri?) and l/{roR). 

C. Case 3. Large Impact Parameter 

1. Asymptotic expansions of independent variables 

In this limiting case we assume that the distance R between observer and source of light is much smaller than r 
and ro, their respective distances from the deflector (see Figure 0). Then we have 



2R r^ 
rl=r'' - 2rRcose + R'' = ( 1 cos6l + — 



which leads to the expansions 



(196) 



ro — r ~ R cos 9 + . 

1 I R 

— = - IH cos ( 

ro r \ r 

The time parameters are 

T ~ r cos 6 , and tq ^ t — R . 
Their numerical values may be larger or less than zero. The following exact equalities hold: 



(197) 
(198) 

(199) 
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{yr) ^ 

In addition, we have asymptotic expansions 



d = r sin , 

y = T — r — r(cos — 1) , 
1 



r'^{cos9 — 1) 



2/0 = To - ro = y ( 1 + — 



(2/0 r-o) ^ 



1 R 

— + — 

yr r"^ 



to - ro = t - r + i?(cos 6-1) 



(200) 
(201) 

(202) 



(203) 

(204) 

(205) 
(206) 



Thus, we can decompose any function of the time argument to — ro in a Taylor series with respect to the retarded 
time i — r if convergence is assumed [ |62[ . For example, 

li] {to - ro) = lij {t-r)+ R{cos 9-1) (i - r) + ... . (207) 

Finally, we note that the vector ^* corresponding to impact parameter d can be represented as 

C ^r{N' ~k' cos 0) , (208) 

where iV* = —Kq = x^/r, |N| — 1, and /c* is the unit vector in the direction from the source of light to observer |]63|] . 

2. Asymptotic expressions for time delay and the angle of light deflection 

In this section all asymptotic expressions for relativistic effects are given only up to leading terms of order 1 /r and 
l/ro- For this reason all residual terms of order 1/r^ and 1/rg are omitted in subsequent formulae without note. 
Using asymptotic expansions of functions from the previous section and reducing similar terms we obtain 



1 



1 — cos ( 



1 



- 2k' cos e + - (1 + cos^ 0) N'N^ 



iij{t-r) jij(io-ro) 



ro 



(209) 



where cos = k • N = K • Ko (see Figures |^ and |4|) . We note that the expression for time delay given above can be 
further simplified if the definition of "transverse-traceless" tensor with respect to the direction A^* is applied (24| , Q : 



Ji^ = J,, + - ((5y + N,N.j) NpN, Jp, - [S^pNjNg + 5,pN,Ng) I, 



PI 1 



(210) 



where the projection is onto the plane orthogonal to unit vector iV\ Formula ( |209D for time delay now assumes the 
form 



1 — cos ( 



if;^{t-r) i,7(to-ro) 



(211) 



Differentiation of Ag with respect to time gives the frequency shift due to a remote localized source of gravitational 
waves 



dt 



(212) 



where the metric hfj^ is defined by the equation (BIC) and taken in the leading order approximation with respect to 
l/r. We recognize that the expression (212) is a generalization of the analogous formula for Zg obtained previously 
by Mashhoon & Seitz |Q in the case of a plane gravitational wave. This exact coincidence demonstrates the power 
of our formalism, which both reproduces well-known results and yields new observational predictions for relativistic 
effects in the propagation of light rays in the field of an arbitrary source of gravitational waves |63] . 
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Repeating the calculations for the angle of light deflection under the assumption that the wavelength, A, of grav- 
itational waves emitted by the localized source is smaller than the distance R between source of light and observer, 
we come to the following result: 



ar 



— - [(cos6i- 2) (k'kPk'' + 2k' cos 6) + (cos^ 6* - 2 cos6i - l) x (213) 

1 — cos '- ^ ' ^ ' 



X ( i/c^iVPiV'cose - N'NPN"^ ] + N'^k^k'^ - iN^N^k'^ 



+2 {kP-NP cose) 



%p[t - r) 



Transformation of this result using relationship ( |210| ) and expression (BIC) for hj-'" ^ where only leading terms of order 
1 /r are retained, reveals that 

1 kPk'i 

= 2 l-k-N ■ ^ " + ^'^^^^ " + ^'^'^P^^ " ' ^^^^^ 
and, because the vector /3' is small, 

s'^K' + o^Q+-f' + .... (215) 

where the ellipsis designates unimportant terms of higher order with respect to 1/r [6^ ], and we have neglected 
the constant deflection caused by mass-monopole and spin-dipole dependent terms. One sees again that only the 
transverse-traceless component ft,^-^ of the metric tensor appears in the final expression. 



It is worthwhile to stress that the observed optical direction to the source of light given by the formula (215) 



coincides with that which can be obtained by means of VLBI observations. Indeed, it is easy to confirm that equation 



(215) can be re-written as follows 



+ \ ^^^^ KPK^hl^it -r)- \KPhl^{t - r) . (216) 



The direction to the source of electromagnetic waves measured by VLBI is determined as difference between times 



of arrival of the wave to the first and second antennas. Taking into account equations (^4|) and (211) for the first 
and second observing sites, and assuming that the time difference t2 — ti in observation of the radio signal at the 
observatories is small compared to the period of gravitational waves, we find 

h~h = -iK + - KPK^hl^it - r)\ ■ (X2 - xi) . (217) 



If the baseline vector measured in the local inertial frame is denoted as b and the transformation (|106|) is taken into 
account, 

x\ -x\=b'-^ h^J{i - r)V + 0(b2) . (218) 

We confirm that 

t2-ii=-s-b, (219) 
where the vector is given by formula (^16|), which proves our statement. It is worth emphasizing that equation 



(216) was obtained independently by Pyne ei al. ( ]16|, see formula (47)). Their formalism, however, has a limited 



region of application. Extension of the formalism of Pyne et al. [|16[ was one of the motivation of the present work. 
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VIII. CONCLUSIONS 



The most accurate astrometric measurements are differential. They measure the angle between 2 sources. The 
highest accuracy is attainable when the sources are close to each other in the sky. In contrast, angular deflection by 
gravitational waves varies only over large angles in the general case of large impact parameter. Specifically, in such 
a case the bendi ng a ng le de pends only on the metric in the neighborhood of the observer and its first derivatives, 
as in equations ( 214 ), ( 215 ). It thu s can vary only as a quadrupole and the derivative of a quadrupole, over the 
sky. Similarly, equations (123), (16C) and ( |195| ) depend on th e mass quadrupole moment lij and its first and second 
derivatives. Note that the the angle of light deflection (123) involves the time integrals oiXijit — r) which may be 
interpreted as the presence of the "kinematic resonance effect" iQ; however, this term is small, as discussed above. 
In the context of a purely locally-determined deflection angle, it is not unexpected that lines of sight that pass close 
to the deflector show almost purely the static effect, as was shown in section 7. 

The magnitudes of the leading terms in the limiting forms for the deflection angle ag, in equations (123), (16C) 
and ( |195D are aq ~ Vt^GMa? /c^r, where M is the mass of the deflector, and a is its dimension. The frequency of 
the gravitational waves is 17. For a gravitationally bound binary system with a circular orbit, f2 is twice the orbital 
frequency |6^. We can use Kepler's third law to express this in the form aq ~ ^^/^{GMY'^^ /c'^r, or alternatively, 

aq ~ 2.4 x lO^^'^(M/M0)^/^Psec^'^(?'kpc)~^ arcsec where Psec is the orbital period of the binary system. For a contact 
white-dwarf binary at 200 pc, the expected deflection is about 7 x 10"^^ arcsec, with a period of about 1000 sec. For 
a supermassive black-hole binary, with mass 10^ M© and period 10 yr at a distance of 1 Mpc, the expected deflection 
is about 5 x 10~^^ arcsec. 

Because the effect varies smoothly over the sky, the presently available astrometric accuracies are a few microarcsec- 
onds. Higher accuracies are attainable only over smaller angles. Very-long baseline interferometry of a suite of radio 
sources attains microarcsecond accuracy, over periods of days to years. Specially-designed observations sensitive to 
source motions of minutes or hours might attain higher accuracy, perhaps as much as an order of magnitude better. 
Clearly, detection of deflection of light rays by gravitational waves from nearby localized sources is not a goal for the 
near future because of its smallness. However the background gravitational wave noise may be, perhaps, measured. 

The near-perfect cancellation of the effect in General Relativity suggests that deflection of light by gravitational 
waves could be a test of that theory in radiative regime |7^ . In a theory that does not posess the symmetries that 
cause the deflection to vanish, we can only guess the resulting deflection. Such a guess might multiply the general- 
relativistic Oq by 3 factors. The first factor, of r/d, reflects the amplitude of the gravitational wave at the point 
of closest approach, rather than at the observer. The second factor, some function of the distance to the source 
measured in gravitational-wave wavelengths, perhaps ln(r/A) ||7l|, reflects the cumulative effect of bending along the 
line of sight. The final factor, unknown, reflects the coupling of the non-general-relativistic part of the wave to the 
source and its effect on the light ray. For a source an arcsecond from the deflectors described above, the first 2 
factors can increase the effect by several orders of magnitude. Moreover, if the effect is not local to the observer, 
differential astrometry across small angles can detect it, so that greater accuracy is attainable. Given sufficiently 
strong departures from General Relativity, the effect might be detectable. 
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APPENDIX A: COMPARISON TO THE PAPER BY DAMOUR &i ESPOSITO-FARESE 



In this Appendix we rederive the results of the paper by Damour & Esposito-Farese |23| applying the generalized 
isotropic ADM coordinate conditions (For an application of the conditions in post-Newtonian calculations, see e.g. 
PH). This explicitly shows that the asymptotic results do not depend on the chosen gauge. We do not put c = 1 in 
this and the following appendices to make more clear the order of terms with respect to the small parameter 1 / c. 

The ADM coordinate conditions, in linear approximation, read 
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2Vi5oi - Voffij = , 3Vj.gjj - Vigjj = , (Al) 

where Vq — d/dt and — d/dx^. For comparison, the harmonic coordinate conditions, in hncar approximation, 
read: 

2Vi5ioi - Vo5ii = Vo5oo , 2Vj.g.y - Vj.gjj = -V^goo ■ (A2) 



The ADM gauge conditions (Al) brings the space-space component of metric to the form 

1 



9^o ^ + ^hkk) + hf^ , (A3) 

where hj^ denotes the transverse-traceless part of hij. Furthermore, in Hnear approximation, the Einstein field 
equations read 

/ioo = -^A-i(Too + r„;) , (A4) 

= -^A-i (^To, - ivoV.A-iToo) , (A5) 

/ifcfc =-^A-iroo , (A6) 

hl^ = -^^P.,Mn^},Tu . (A7) 

where the TT-projection operator Pijki, applied to symmetric tensors depending on both time and spatial coordinates, 
is given by 

P^m = {5^k - A-iV,;Vfc)(<5,, - A'^V i) - ^(^ - A'^V ,V ,){5ki - A-iVfcV,) , (A8) 

where denotes the Euclidean inverse Laplacian. 

We now follow the calculation in the paper by Damour & Esposito-Farese |Q using the ADM coordinate conditions. 
The deflection of the light ray is given by 

= - l^i-^) = - / drl^fy ^^KpiS,^ + Tl^). (A9) 

^ J -oo 

In terms of the Fourier transform 

V(fc^)= J d^xh^.(x^) e~''''-\ (AlO) 
where fc^x" = —k^x^ + k • x, the boldface lett ers d enote spatial components of vectors [Q, and integration is over 



all of space-time. Accounting for the formulae (AlO) the equation for light deflection now reads 

Al^=in J e^"-«"r/^L/3(A:^)<5(k.l-fc°Z°), (All) 

where use has been made of the exact relationship 

/ + 00 
cLt expiirkj") = 27r5 (k • 1 - k"l") . (A12) 
-oo 

In terms of the Fourier transformed energy-momentum tensor Tap, the Fourier transformed metric field reads 
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(A13) 
(A14) 
(A15) 
(A16) 



where e is a positive infinitesimal number which shows exphcitly that we have used the retarded Green's function 
while solving the Einstein equations for the component of the metric hij. If we put k • 1 — k^P = and take into 
account the relations 



rp I- J rp 



Tm — -r-Ti 



ko 



which follow from the macroscopic equations of motion for matter V i,T'^'^ — 0, we find: 



h 1 



T- 



(A17) 



(A18) 



This expression is identical with that obtained by Damour & Esposito-Farese in the harmonic gauge. For this 
reason calculation of the total deflection angle gives the same result in both harmonic and ADM gauges, reflecting 
the coordinate independence of the final result. 

Defining = Al^^/l'^, one gets for the angle of total deflection 



rf3 



"2 = -^^12{t*), 



a,=a =-- 



iii(r)-i22(<*) 



(A19) 
(A20) 
(A21) 



One can see that expressions for ai, a2 are the same as those obtained in section 7.1.2. The quantity a*^ gives the 
gravitational shift in frequency of the electromagnetic wave. It can be obtained from the expression for gravitational 
time delay Aq after its differentiation with respect to time. 



APPENDIX B: HARMONIC AND ADM GAUGE CONDITIONS IN THE FIRST POST-MINKOWSKIAN 

APPROXIMATION 



In this Appendix we give other representations of the metric coefficients ( |12D - . Using the ADM coordinate 
conditions ( |Al| ) of Appendix A the metric coefficients (12) - (|l^ can be cast into the "canonical" ADM form 



^So^'" = ^ + ^V.V,.-\ (Bl) 



"-Oi - ^3 ^2 + ^3 4^3 VjVjVfcr, [HZ) 

h-kkT = 3 /loo, (B3) 
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(B4) 



Iki{t-r/c) Iki{t) 



r r 2c? 

where we have used a special symboHc notation for "semihereditary functionals" | 



(B5) 



The foUowing equahty holds: '^'^iij{t — r) = Tij{t — r). We also notice that A{lij{t — r/c)/r) — iij{t — r/c)/c^r 
for r ^ 0. For this reason function '-^^^Xij{t — r) is a solution of the homogeneous d'Alembert's equation; that is, 
□ ['^-^^I^J{t - r)/r] = for r ^ 0. 



We emphasize that the metric (B1)-(B4) is an external solution of the equations (A4)-(A7), outside the source of 



gravitational waves. It matches smoothly to the internal solution, which is valid inside the source, without additional 
coordinate transformations. It is remarkable that outside the source the metric component (B4) may be represented 
as an algebraic decomposition of the retarded and instantaneous functions of time 



h: 



admTT 



(B6) 



where hjj'"{t — r/c, x) is shown below in ( BlOj ), and hij{t,x.) is the rest of the metric component /ifj*™-^-^ which 
is actually a symmetrized gradient of a vector comprising of singular harmonic functions. For this reason function 
hij(t, x) satisfies the condition Piju hki{t, x) = and can be eliminated by an infinitesimal coordinate transformation. 
Making use of this and without leaving the ADM coordinate conditions, we may construct the following representation 
for the metric: 



/loo 
hkk 
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(B7) 

(B8) 
(B9) 
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1ki{t - r/c) 



^-^)lki{t-r/c) 



This form of the metric is obtained from the expressions (B1)-(B4) by applying the coordinate transformation fl_ 



w 



cr 



2Vfc 



ikiity 



(Bll) 



(B12) 



It is marvelous that this representation of metric also fulfils the harmonic coordinate conditions (A2). This means 
that outside the localized source of gravitational waves the class of ADM coordinates overlaps with that of harmonic 
ones. The coordinate t rans fo rmat ion from metric (12)-(^4|) written in the pure harmonic coordinate system to the 
ADM-harmonic metric (^)-( |B10| ) is: 



(-i)X,i(t-r/c) 



1. 



w 



'-y^ykyl 



cr 

(-2)Xfe,(i-r/c) 



2Vfc 



^ki{t - r/c) 



(B13) 
(B14) 
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These gauge functions have been extensively used in the main body of the paper for elaborating unique interpretation 
of observable effects. In contrast to the expressions (B1)-(B4) the expressions (B6)-(B9) show terms which decay like 
l/r* and l/r^. These terms depend on time integrals of the quadrupole moment lu (t — r/c) and may lead to the 
appearance of the "kinemati c res on ance effect" discussed by Braginsky & Grishchuk [Q . Another important remark 
is that the transformations (B13)-(B14) clearly show how to eliminate non-radiative terms from the metric, written 
down in a harmonic gauge, including all terms with respect to any power of 1/r. Previously used transformations 
(see, for example, the textbook of Misner et al. [Q, paragraph 35) dealt only with terms of the first order in 1/r 
and could not be applied for analysis of gravitational radiation in near or intermediate zones of the localized source 
of gravitational waves. 
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itational wave detectors such as LIGO or LISA. For instance, if we introd uce the notation AR = cAq for the variation 
of distan ce b etween the test bodies of the detectors, then from formula (^11) and Taylor expansion of the quadrupole 
moment (207) one immediately obtains the well known result 
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where A is the wavelength of gravitational waves impinging on the detector. If necessary we can write down corrections to 
this formula up to any order. The case when A is comparable with R should be treated in a slightly different way and is a 
matter for a separate study. We emphasize that this can be done in the framework of our formalism as well. 
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[69] Tsvetkov, V. P., & Tsyrulev, A. N., 1983, Sov. Astron., 27, 66 

[70] Damour, T., & Esposito-Farese, G. 1992, Class. & Quant. Grav., 9, 2093 

[71] Faraoni, V., & Gunzig, E., 1998, A&A, 332, 1154 

[72] For example, k = it* = (fc\ fc^ fc^). 

[73] Blanchet, L., & Damour, T. 1988, Phys. Rev. D , 37, 1410 

[74] We recall that /ifi^ h^v + V^w^ + V^Wi^ and indices are raised and lowered with the help of the Minkowski metric rjap = 
diag(— 1, +1, +1, +1). The spatial part of the transformation Vi^iu^ + Vpiu,/ also satisfies the identity Pijki ( VfeiOj + VjWfe) = 
as does the function hij{t,x). F or thi s reason the TT-part of the metric (B4) can be transformed into (BIO) by means 
of the coordinate transformation (B12). This explicitly reveals that the decomposition of the spatial component of the 
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FIG. 1. Astronomical coordinate system used for calculations. The origin of the coordinate system is at the center-of-mass 
of the source of gravitational waves. The bundle of light rays is defined by the vector field fe*. The vector = — fc' + 0(c~^) 
is directed from observer towards the source of light. The vector Kq is directed from the observer towards the source of 
gravitational waves. Wc use in the paper the equalities Kl) = —A''' = jr, where arc the coordinates of the observer with 
respect to the source of gravitational waves, and r = |x|. The plane of the sky to the vector K}^ is not shown. 
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FIG. 2. Relative configuration of observer (O), source of liglit (S). and a localized source of gravitational waves (D). The 
source of gravitational waves deflects light rays which are emitted at the nioincnt to at the point S and received at the nionient 
t at the point O. The point E on the line OS corresponds to the moment of the closest approach of light ray to the deflector 
D. Distances are OS = R, DO = r, DS = ro, the impact parameter DE = d, OE = r > 0, ES = tq = t — R < 0. The impact 
parameter d is small in comparison to all other distances. 
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FIG. 3. Relative configuration of observer (O), source of light (S), and a localized source of gravitational waves (D). The 
source of gravitational waves deflects light rays which are emitted at the moment to at the point S and received at the moment 
t at the point O. The point E on the line OS indicates the point of minimal distance of the light ray trajectory extrapolated 
backward to the deflector D. Distances are OS = R, DO = r, DS = ro, the impact parameter DE = d, OE = r > 0, 
ES = To = T — R > 0. The impact parameter d is small in comparision to all other distances. 
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FIG. 4. Relative configuration of observer (O), source of light (S), and a localized source of gravitational waves (D). The 
source of gravitational waves deflects light rays which are emitted at the moment to at the point S and received at the moment 
t at the point O. The point E on the line OS corresponds to the moment of the closest approach of light ray to the deflector D. 
Distances are OS = R, DO = r, DS = ro, the impact parameter DE = d, OE = t = rcoaO, ES = to = r — R. The distance 
R is much smaller than both r and ro- The impact parameter d is, in general, not small in comparision to all other distances. 
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